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Abstract Surface waves on both superfluid 3 He and 4 He were examined with 
the premise, that these inviscid media would represent ideal realizations for this 
fluid dynamics problem. The work on 3 He is one of the first of its kind, but on 
4 He it was possible to produce much more complete set of data for meaningful 
comparison with theoretical models. Most measurements were performed at 
the zero temperature limit, meaning T < 100 mK for 4 He and T ~ 100 K 
for 3 He. Dozens of surface wave resonances, including up to 11 overtones, were 
observed and monitored as the liquid depth in the cell was varied. Despite of 
the wealth of data, perfect agreement with the constructed theoretical models 
could not be achieved. 

Keywords Surface wave • Superfluid • Normal fluid • Interdigital capacitor 


1 Introduction 

Fluid dynamics is an important branch of physics. However, accurate modeling 
of fluid flow in general geometries is cumbersome, and proper account for the 
effects of viscosity often leads to serious difficulties. Superfluids, on the other 
hand, offer the possibility to study quite common flow problems without the 
complications due to viscosity and, thus, they represent the closest possible 
realization of the ideal fluid often hypothesized in theoretical models. Also, 
the extraordinary wetting properties of the superfluids lead to well defined 
conditions in cases, where the free fluid surface comes to contact with solid 
walls. For these reasons, it was thought that the surface waves on the free 
surface of superfluids would obey a simple theoretical description nice and 
easy. In practice, however, it turned out not so straightforward to obtain good 
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agreement with the experiment and the theory. This paper summarizes our 
understanding of the problem at present. 

Surface waves on free fluid surfaces are perhaps the most easily perceived 
representation of wave phenomena in our common environment. The standing 
wave modes in restricted geometries can easily be demonstrated and analyti¬ 
cal solutions to these in some idealized cases have been known for long. The 
resonance frequencies of the surface modes are influenced by kinetic, potential, 
and surface energies, and the problem in question can be viewed from several 
perspectives. Both bulk fluid and surface properties are connected together 
and they need to be incorporated to any proper model. It must be acknowl¬ 
edged that the free fluid surface in a bound geometry is not perfectly flat, but 
the menisci caused by the surface tension have an important contribution to 
the energy balance. Under general circumstances the vapor phase above the 
fluid surface would have a minor contribution to the problem but in the case 
of helium fluids at very low temperatures the gas phase becomes extremely 
rare and can be ignored completely. 

Even if we are dealing with superfluids here and the fluid itself can flow 
without friction at low velocities, the surface waves still experience damping 
at any finite temperature. This is due to interaction of the moving interface 
with thermal quasiparticles in the fluid, and can be measured by observing 
the frequency width of the surface wave resonances. At sufficiently low tem¬ 
peratures the quasiparticles become ballistic, they scatter from the surfaces 
only, and their effect can be formulated easily. Interestingly, the quasiparticles 
have quite different character in bosonic 4 He and in fermionic 3 He superfiuids, 
which lead to entirely different temperature dependencies. This feature has 
been presented in more detail in an earlier paper T], but will be reviewed 
shortly here again. 

We studied standing waves on free fluid surface both in superfluid 3 He and 
4 He. Surface waves in superfluid 3 He have been studied only very recently [2] 
thus our work in 3 He pushes to new frontiers with this rather exotic substance 
and the required microkelvin range temperatures. However, our data on 4 He 
is much more systematic and complete, and thus better suited for meaningful 
comparison with theoretical models. Nevertheless, interesting comparison can 
be made between the results on those two fluids investigated in the exact same 
experimental cell. 

In this paper our focus is mainly on the resonance frequencies and energy 
balance in the zero temperature limit, though we give some remarks on the 
temperature dependencies of the characteristic quantities, and shortly discuss 
the surface waves above the superfluid transition temperature T c . We were able 
to distinguish and identify several resonance frequency modes with varying 
depth of liquid helium in the system. The resonance frequencies were deter¬ 
mined with great precision but the nontrivial geometry of our cell prevented 
us from deriving a conclusive value for surface tension in the zero temperature 
limit. Selected frequency bands could be fitted with credible values from the 
literature but the best overall correspondence across the whole data set led to 
somewhat contradicting values for some system parameters. 
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Fig. 1 (Color online) Geometry of the experimental cell in top view (left) and side view 
(center). The annealed copper cell walls (dashed) were coated with about 0.5 mm thick 
layer of sintered silver powder heat exchanger (gray). The two interdigital capacitors (IDC) 
were located at the opposite vertical walls of the cuboid volume. The full sensitive range 
of the IDCs is 4.6 mm in vertical direction and it starts at about 1.4 mm from the porous 
sintered bottom (or 2.0 mm from the impermeable copper bottom, where the thickness of 
the sintered layer is about 0.6 mm). The two capacitors were slightly tilted, as indicated by 
the measured angles. The IDCs were thermalized to the cell wall with copper mesh (right) 


2 Experimental setup 

Oscillations of the free surface of liquid helium were studied in an experimen¬ 
tal cell illustrated in Fig. [TJ Heat was transferred from the helium sample via 
sintered silver powder heat exchangers to the copper nuclear cooling stage of a 
microkelvin refrigerator [3|. Temperature was measured at different stages by 
a germanium resistor, a SQUID-based noise thermometer, a platinum NMR 
thermometer, or deduced from the demagnetization field at the lowest tem¬ 
peratures. 

The fluid surface level was sensed with an interdigital capacitor (IDC) on 
a vertical wall of the cuboid volume with time resolution of about 10 ms, al¬ 
lowing detection of resonances below about 100 Hz frequency. The detection 
signal amplitude had to be kept low for two reasons: first because large readout 
signals generated heat, which was detrimental for the experiment, in particular 
in the case of superfluid 3 He, but also because large electric field/heated ca¬ 
pacitor surface attracted superfluid on the surface and altered the conditions 
of the measurement. Level sensitivity of about 0.1 /im was achieved without 
notable undesired effects. Another capacitor was located on the opposite wall 
of the cuboid by which the surface waves were meant to be excited. How¬ 
ever, excitation by means of electric field was not successful and eventually 
the waves were excited mechanically either by ambient vibrational noise or by 
periodically swinging the whole cryostat in micrometer scale. The second IDC 
was then used to monitor the static fluid level in the cell. The wave resonances 
generated by ambient vibrational noise were found from the Fourier transform 
of the IDC data processed by an oscilloscope. Alternatively, when the waves 
were generated by swinging the cryostat with a series of different frequencies, 
the IDC signal with 10 ms time resolution was further processed with a lock-in 
amplifier. For details, see Ref. [4j. 
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The IDCs were microfabricated elsewhere by A. J. Niskanen, Aalto Univer¬ 
sity, Microfabrication group. Sapphire substrate (0.3 mm x 10 mm x 10 mm) 
cut orthogonal to the c-axis was patterned with 740 micro strips (spacing 
5 fj, m, width 5 /im) and the effective area of the IDC was 7.4 mm wide and 
4.6 mm high. The patterned deposited layer of aluminum was 50 nm thick with 
additional 15 nm layer of chromium to improve adhesion. The pattern was su¬ 
perconducting below 1.3 K. Each IDC was thermalized with a copper mesh 
(see Fig. |I| which was glued (Stycast 1266) to the back of the sapphire plate 
with an extension bolted to the cell wall. Superconducting niobium-titanium 
wires (filament diameter, 62 fim; bare diameter with copper matrix, 100 fxm.) 
were glued to the copper mesh, and the superconducting filaments whittled 
out under the copper matrix were bonded to the contact pads of the IDC with 
multiple aluminum wires. The thermalization mesh passed through the narrow 
channel from the cuboid volume to the annular volume making the geometry 
somewhat ill-defined there. 

The experimental setup was originally designed for studies at the melting 
pressure of helium [5]. Two quartz tuning fork resonators, usable as thermome¬ 
ters in liquid helium, were thus located at the upper part of the experimental 
cell, whereas the lower part of the cuboid volume with the IDCs was origi¬ 
nally reserved for solid helium. The helium crystal size could be varied with 
a separate bellows volume with base temperature not lower than 0.4 mK. A 
pressure gauge was also attached to the bellows volume. The experimental 
cell and the bellows volume were connected together with a tube (diameter, 
2.5 mm; length, 20 cm). The bellows volume and the quartz tuning forks are 
not shown in Fig. [I] as they were not in use during the experiments with free 
surfaces. However, if positioned differently, both the bellows volume and the 
quartz tuning forks would have been useful assets during the experiments with 
free surfaces also. 


3 Surface Wave Resonances in Superfluid 4 He 

We first present here the observed surface wave resonance modes in superfluid 

4 He at the zero temperature limit (< 100 mK) as the function of liquid depth 
in our experimental cell, as driven by the ambient vibrational noise. This is our 
most complete set of data, which we then analyze in terms of the theoretical 
models in the following sections. Auxiliary data on the behavior as the function 
of temperature, the data on superfluid 3 4 He, and some observations in the 
normal fluid He are presented in the latter sections. 

An example of the Fourier transformed signals is presented in Fig. [2] Every 
such trace is an average of about 20 recordings, each collected for 50 s. The 
noise bottom had slight dependence on frequency but it remained pretty well 
constant for the whole duration of the experiment, which took about 24 days. 
Above the background in the example spectrum, one can observe some sharp 
peaks below about 10 Hz, at about 32 Hz, 40 Hz, and at about 60 Hz, the 
50 Hz contamination at the line frequency and a broader feature slightly above 
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FFT (V 2 /Hz) hi DC (mm) /i IDC (mm) 


Fig. 2 Surface wave resonances in superfluid 4 He generated by ambient vibrations in 
the zero temperature limit (T < 100 mK). Left: A raw FFT spectrum at helium level 
/iidc = 4.4 mm with fitted background. 1 V corresponds to 0.17 mm on the IDC. Center: 
Background subtracted FFT spectra with grayscale encoded intensity at various helium lev¬ 
els hiDC- Right: FFT data after filtering and image processing. Only the response due to 
the liquid surface modes remains visible. The resonance modes corresponding to the cuboid 
experimental volume are labeled with indices ( m x ,m y ) 
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55 Hz. From an individual sample it is not possible to tell, which of those arise 
from the oscillations of the helium surface, and which were due to electric 
disturbances or other artifacts of mechanical origin. It is the continuous filling 
of the cell and the consistent shift of the relevant resonances as the function of 
the liquid level, which tells us the difference. We can thus further clean up the 
data by eliminating any features above the background, that do not depend on 
the liquid level in the cell. The middle panel in Fig. [2] shows the background 
subtracted signals as the function of the liquid depth, and the right panel is a 
cleaned version, where only the interesting features due to the modes on the 
liquid surface remain. 

Such processing reveals a lot of details not distinguishable at individual 
spectra. There are about a dozen surface resonances below 15 Hz frequency 
and a more or less evenly spaced set of 12 resonances up to 60 Hz frequency. No 
further modes could be observed above that. It is noteworthy that the liquid 
level dependency of the modes at the highest frequencies was almost entirely 
due to the slight tilting of the capacitor plates, since the surface dimensions 
then slightly depended on the liquid level. Were they perfectly aligned, it 
might have been impossible to distinguish these modes from the background 
disturbances! The analysis of these sets of resonances is presented after the 
introduction of the theoretical models in the next section. 

The broad feature in the example spectrum in Fig. [2] at about 55 Hz may 
deserve a special note. This ” bump” was nearly always visible in the data but 
its position varied seemingly arbitrarily in the course of time from about 20 to 
55 Hz. This caused the wandering shadow, that is visible in the middle panel 
of Fig. [2j most clearly at about 3 mm level below 45 Hz and above 4 mm 
level below 55 Hz. We have no explanation on its origin but it was definitely 
not due to any action of the superfluid in the experimental cell. The vertical 
shadows in the middle panel of Fig. [2] were caused by slight variations of the 
noise bottom of the spectra. It also appears that the ambient vibration levels 
at our experimental location were exceptionally high during the measurements 
at about 1 mm liquid level, as the signals from the motion of the surface appear 
more intense as usual during that period of time. 


4 Models for Surface Resonances 

Standing waves on the free fluid surface are modeled assuming incompressible, 
irrotational, homogeneous and inviscid fluid. The parameters characterizing 
the geometry of the cuboid volume are shown in Fig. [3] In Cartesian co¬ 
ordinates the free surface profile can be written in the form z SUI f(x,y,t) = 
h + r](x,y) + £(x,y,t), where y(x,y) describes the stationary meniscus and 
£(x,y,t) = ((x,y)r(t ) represents a standing wave mode on the surface with 
r(t) = asin(o;t). For a semi-infinite horizontal fluid surface bound by a wall, 
that is aligned with the y-axis, the shape of the meniscus r][x) is obtained 
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4.6 mm 


^IDC 


Fig. 3 Free surface in the cuboid vol¬ 
ume. The lowest measurable helium level 
is h = ho when time — 0 mm. For the 
highest measurable level /iidc — 4.6 mm. 
The capillary raise is y/2l c = y/2^ / {pg). 
Volume width at the bottom is L x q 


from [ 6 ] 



(i) 


where l c = y/ 7 / (pg) is the capillary length, and 7 , p and g are the surface 
tension, mass density and gravitational acceleration, respectively. The param¬ 
eter a"o is included for setting the proper contact angle with the fluid. For a 
vertical wall with zero contact angle f r/dx = 1%. 

For the surface tension in 4 He we use 7 = 375 pJ/m 2 [7] and thus l c = 
0.513 mm. For 3 He we use 7 = 155.7 pJ/m 2 [18] and l c = 0.440 mm. 


4.1 Basic Analytical Model 


Let us begin with a basic model, where the meniscus effect on the boundaries of 
the free surface is neglected. Under that assumption the low amplitude stand¬ 
ing surface wave resonances in a rectangular liquid pool occur at frequencies 
/ = uj/(2tt) as 


2 

-°'basic 


k 2 



gk tanh(fcft), 



m x ,m y £ No, 


( 2 ) 

( 3 ) 


where h, L x , and L y are the liquid depth, and the surface length and width, re¬ 
spectively. Further, we have the surface wave profile z sur {(x, y, t) = h+£(x, y , t), 
the scalar potential <f>(x,y,z,t) [defining the velocity field v(x,y,z,t) = V0], 
and kinetic and potential energies -Ekin and E pot as follows: 


£(x,y,t) = cos 


(f>(x, y, z, t) = cos 


-Ukin — 


m x ir 


x cos 


TO r 7T 


X COS 


V T (t) = COUZ/M t), 


( 4 ) 


J v 


cosh (kz) 1 
sinh (kh) k 


■(t) = ^(x,y,z)f(t),( 5) 


1 2 7 T L X LyP 

— o\ dv =--- 

2 P 2 ktanh(kh)2 D 


= ~M kin , 


( 6 ) 
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-E'pot — 


;P<A 2 


(W + (dyZY 


dx dy 


_ t 2 uj 2 L x L y (pg + yfc 2 ) _ t 2 lo 2 

2 u 2 2 D 2 pot ’ 


(7) 


where the time dependencies are explicitly separated as the factor r(t) = 
asin(wf). The parameter D depends on the indexes m x and m y : having either 
m x = 0 or m y = 0 implies D = 1, and with both m x > 1 and m y > 1 one has 
D = 2. From the Euler-Lagrange equation, r being the generalized coordinate, 




dL 
dr ’ 


L — Ek in — E , 


pot 


( 8 ) 


it follows that M = A/kin = M pot for the harmonic oscillator -type equations 
Eqs. ® and 0, in agreement with Eq. ©■ 

Similar equations can also be written for the annular volume in cylindrical 
coordinates. However, it is simpler to describe such a narrow annular volume 
with the equations derived for the rectangular geometry without any loss of 
accuracy — only the boundary condition in one direction becomes periodic 
due to the solutions wrapping up to themselves around the circumference. 


4.2 Analytical Model with Meniscus Correction 


In order to account for the meniscus effect Roche et al. [7] treated the meniscus 
r)(x, y) in a wide well ( m y = 0, L y L x ) as additional mass on the free surface 
with the profile x SU rf(x, y,t) = h + r](x , y) + £(x, y)r(t). Surface displacement, 
scalar potential and thus also the total potential energy are assumed to be of 
the same form as Eqs. 0 ,® and 0 without the meniscus corrections. These 
assumptions are valid only when kl c 1 : the meniscus height i) ~ Z c should 
be smaller than the depth of the layer where the velocity field is significant, 
which is given by A: -1 . However, kinetic energy increases due to the moving 
additional mass in the meniscus as 


Ekin 


t L x L y p 
2 2 fctanh(fc/i) 


<7 1 

1 + 4 k -tanh (kh) 

P9 Lx 


(9) 


and the Euler-Lagrange equation Eq. (| 8 ]) can be used to obtain the resonance 
frequencies 


UJ 


2 


_^basic_ 

1 + 4fc ^^ tanh ( fc/l )' 


( 10 ) 


The original result by Roche et al. dealt with just one of the two menisci on a 
straight channel whereas the above equation takes both of the opposite walls 
into account. 
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Fig. 4 (Color online) (Bottom) 
scalar potential ip(x, z) calcu¬ 
lated with FEM in 2D accord¬ 
ing to Eqs. (113D , (11411 and (1 15 D . 
(Top) the black line illustrates 
the exaggerated free surface de¬ 
viation from the equilibrium in 
the normal direction n. The si¬ 
nusoidal form is shown for com¬ 
parison as the thin gray line. 
The mesh for the model was 
generated with aid of the code 
described in Ref. [9] 


In the same spirit as above it is possible to include also the two additional 
menisci on the side walls in the case when m y = 0 but both L x and L y are 
finite: 


k^basic 


l + 4/c-£ 

P9 


tanh (kh) + [tanh (kh) + coth(/c/i)]| 


( 11 ) 


In an annular volume menisci exist on the sides only; thus for the case with 
periodic boundaries one gets 


UJ 


2 


_^ basic_ 

1 + 2/c^-^ [tanh(/c/i) + coth(/c/i)] 


( 12 ) 


4.3 Numerical Model with Meniscus 

The model described in the previous section is limited to wavelengths much 
longer than the capillary length. To extend the validity of the analysis to 
shorter wavelengths and for comparison with cases solvable by other means we 
have performed numerical computations by the finite-element method (FEM). 
It is then straightforward to deal with the geometry altered by the menisci, but 
it is not so obvious, how to incorporate all effects of the surface tension into 
the system of equations to be solved. We follow a partly phenomenological 
approach, where the modified geometry due to the surface tension is built 
into the model to begin with but the influence of the surface tension to the 
energy balance is incorporated as a scaling factor for the potential energy in 
the system. This results in eigenvalues scaled by factors depending on the full 
solution of the wave profile. This approach was proven valid in cases, that can 
be compared with analytical solutions. 

For the starters the spatial scalar potential ip(x, z) and the resonance fre¬ 
quency wo are computed in a geometry with curved free surface from the 
eigenvalue equation 

V 2 V> = 0, 

n ■ S7ip = 0 on solid walls, 


(13) 

(14) 
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n • Vip = —— , =ip on free surface, (15) 

9 y/1 + {d x rif 

where n is the surface normal vector. An example solution corresponding to 
m x = 3 is shown in Fig. © 

As the final step to obtain the actual resonance frequencies u> from the 
Euler-Lagrange equation Eq. ©, surface tension is included as an extra po¬ 
tential energy term, while assuming that the numerically calculated scalar po¬ 
tential remains unchanged. The resulting resonance frequencies w are attained 
from ip and wq as 


2 2 
w = U) Q 



pg 


If surf [1H -(a*)7 ) 2 ] 3/2 ^ X dy \ 

i/surf d X d'lJ )' 


(16) 


Without the meniscus this modeling principle is exact, as the scalar potential 
in Eq. © does not depend on the surface tension at all. The situation is 
not quite as clear, when dealing with the curved surfaces, but the anticipated 
magnitude of the resulting error is not significant. This model will be discussed 
in more detail in the Ph.D. thesis by Matti S. Manninen. 


5 Comparison of Experiments and Theory 

Let us now further examine the observed surface wave modes in superfluid 
4 He studied below 100 mK while the cell was continuously filled up. 

The resonance frequencies corresponding to the geometry of the cuboid 
volume are presented in Fig. [5] Altogether twelve resonances can be identified 
with mode indices from (m, = 1, m y = 0) to (m x = 12, m y = 0). Due to the 
symmetries of our geometry the capacitors were insensitive to the modes with 
m y > 0. The lowest measurable helium level ho = 1.3 mm was treated as a 
fitting parameter. Also the angle 9 = 2.75° between the two capacitors was 
determined from the data. 

Evidently the basic model, Eq. ©, where menisci have been neglected, is 
not sufficient, see circles in Figs. [5] and [6] Analytical 2D meniscus correction 
[Eq. (1101) ] is not any better as it overestimates the kinetic energy related to 
the menisci especially when the condition kl c <C 1 is not fulfilled (triangles in 
Fig.©. Moderate agreement with measurements is reached with a numerical 
2D model, Eq. m , where the meniscus effect is taken into account on the two 
IDC end walls but neglected on the silver sintered side walls (plus-symbols in 
Figs. ©and©. However, the fluid contact angle on the porous sintered side 
walls may also be close to perfectly wetted zero angle and thus the meniscus 
effect on the side walls should not be omitted. To illustrate the effect of the 
menisci on the sides, we have treated them numerically in ^-direction [Eq. (1161) ] 
and in y-direction as an extra kinetic energy term, ^p f 2l^v 2 dx , see cross¬ 
symbols in Fig. © However, this model does not fit the data as the model 
overestimates the kinetic energy in the menisci. Evidently, in order to take 
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Fig. 5 Left: FFT spectra of the surface waves in superfluid 4 He generated by ambient 
vibrations in the zero temperature limit (T < 100 mK) at various helium levels hiDC- Solid 
lines represent eigenfrequencies in the cuboid volume without meniscus correction according 
to Eq. j2j- Dashed lines are based on 2D FEM calculations according to Eq. JTfi}. The 
measured resonance frequencies are highlighted with dots. Frequency differences between 
each model and the data are shown as circles and plus-symbols. Here 7 = 375 pJ/m 2 , 
L x 0 = 10.0 mm, 6 = 2.75°, and Hq = 1.3 mm. Right: Frequencies up to 13 Hz enlarged 
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into account the menisci also on the sides, a full 3D numerical model should 
be introduced. 

Although the basic model, Eq. 0, does not properly fit the data, we can 
estimate the excess energies from the frequency difference between the simple 
model and the measurements. We scale the kinetic energy of the basic model 
Eq. 0 by a factor of (1 + Am x ) whereas the potential energy remains as in 
Eq. 0. As the result, according to the Euler-Lagrange equation Eq. 0, the 
resonance frequencies in Eq. 0 are scaled by the factor of (1 + Ato x ) -1 / 2 . We 
get A = 0.014 for this phenomenological fitting parameter, and the corrected 
resonances are shown as gray filled circles in Fig. [G] The corrected resonances 
are in good agreement with the data for all the twelve modes. 

To justify the chosen form of the scaling factor (1 + Am x ), we note that it is 
of the same form as Eq. ©, triangles in Fig. [0 though A = Att(I c /L x ) 2 = 0.033 
in this case for high helium levels. 

All our models deal with some liquid in a free volume bound by solid walls. 
However, the bottom and the side walls of our cell are coated with a layer of 
porous sintered silver powder. The fitted value for ho = 1.3 mm in Fig. [5] 
suggests that the distance from the effective region of the IDC to the bottom 
is ho + a/2 l c « 2.0 mm, see Fig. [3] The measured value for the distance from 
the effective region of the IDC to the sintered bottom is about 1.4 mm and 
to the impenetrable copper bottom about 2.0 mm. Thus we conclude that the 
sintered walls cannot be considered as solid walls for superfluids. The sintered 
layer is not a free volume either but rather a region with lower mass density 
of the fluid. For the modes with higher frequency or when the helium level 
is sufficiently high, the effect caused by the porous sintered layer should be 
small. Therefore it is not likely that such effects would explain the difficulty 
of fitting our data with the presented models. 

Figs. 0a) and 0b) are shown to demonstrate that our models can repro¬ 
duce the measured data if some model parameters are set free. In Fig. 0a) 
the analytical meniscus corrected 2D model [Eq. (flOl) ] is fitted to the data 



m x 


Fig. 6 Frequency difference between various 
models / mo del and measured resonance fre¬ 
quencies /meas at helium level of /iidc — 
4.6 mm. As in Fig. [5j 7 = 375 pJ/m 2 , L x o = 
10.0 mm, L y = 10.35 mm, and 6 = 2.75°. 
Open circles: the basic model without menis¬ 
cus correction [Eq. J2}]. Gray filled circles: the 
resonance frequencies of the basic model af¬ 
ter scaling by factor of (1 + 0.014m x ) — 1 / 2 . 
Triangles: Analytical 2D meniscus corrected 
model [Eq. dIQD l. Plus-symbols: 2D FEM cal¬ 
culation [Eq. (116ti l . Cross-symbols: 2D FEM 
calculation [Eq. dl6D l with extra mass menis¬ 
cus correction on sides 
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from m x = 2 to m x = 6, when kl c = (m x -K/L x )l c < 1, resulting in the fitting 
parameters L x o = 9.45 mm and ho = 1.22 mm. Another model, where the 
meniscus effect was calculated numerically in ^-direction [Eq. (11611 ] and as an 
extra kinetic energy term in y-direction, is shown in Fig.[7]]b). Also this model 
fits well to the data if we set L x o = 9.0 mm though the measured value is 
about L x o = 10.0 mm. 


5.1 Low Frequency Modes and Temperature Dependencies 

The main cuboid volume of our experimental cell should support only one ob¬ 
servable primary mode of the free surface below 10 Hz frequency. In contrast, 
the low frequency part of our data was decorated by several resonances in that 
region. At least some of those are obviously modes with the largest amplitude 
in the annular volume surrounding the cuboid volume. We can treat the annu¬ 
lar modes according to the argument given in the theory section, although the 
connection with the cuboid volume through the narrow channels reserved for 
the wiring etc. cannot be treated in any easy fashion. Such modes apparently 
produce significant oscillations at the face of the IDCs, though we were not 
able to perform any simple analysis producing good correspondence with the 
actual results in this regard. 

Again, the meniscus effect has to be taken into account. An analytical 
meniscus correction on the side walls is given by Eq. m, where the condition 
kl c <C 1 is better fulfilled for the long annular channel than for the 1 cm 
long cuboid volume. However, in reality the equilibrium free surface profile 
is more complicated and poorly known since the menisci are affected by the 
thermalization grid of the IDCs (two strips of copper mesh passing through 
the annular volume). 

As shown in Fig. El eight annular modes can be distinguished, though the 
overall correspondence is not as good as for the cuboid volume. In addition, 
some modes seem split both in the annular volume and in the cuboid volume, 
see Figs. [8] and El 

We also shortly comment on the artificial periodicity in some of the reso¬ 
nance frequencies, see, as an example, the mode with the second lowest fre¬ 
quency in Fig. El This 0.38 mm periodicity probably reflects the vertical peri¬ 
odicity of the IDC thermalization mesh fed through the annular volume. 

In these measurements we could not observe any temperature dependencies 
of resonance frequencies, which would have reflected changes in surface ten¬ 
sion due to ripplons, quantized surface excitations [TO], see the raw frequency 
spectra in Fig. [9] Instead, our focus was more on the temperature dependency 
of resonance widths Af . Temperature dependencies of Af of three resonance 
modes are shown in Fig. HOlab As expected, the resonance amplitude A is 
inversely proportional to the resonance width Af with constant excitation, 
see Fig. fTOl bb If we assume this proportionality (AAf = constant), then the 
integral of the FFT spectrum over the resonance is inversely proportional to 
the resonance width, J FFT df oc A 2 Af oc ( Af )~ 1 . With the FFT integrals 
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Fig. 7 Measured resonance frequencies in superfluid 4 He (dots). Solid line and circles are 
resonance frequencies according to the basic model in Eq. GQ. (a) Dotted lines and trian¬ 
gles represent analytical meniscus corrected 2D model in Eq. (HOD . Here 7 = 375 pj/m 2 , 
6 = 2.75°, ho = 1.22 mm, and L x q = 9.45 mm. (b) Dashed lines are based on 2D FEM cal¬ 
culations according to Eq. m with meniscus correction in y-direction as additional mass. 
Here 7 = 375 pJ/m 2 , 6 = 2.75°, ho = 1.22 mm, and L y = 10.35 mm, but L x 0 = 9.0 mm 
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Fig. 8 FFT spectra of the surface waves in superfluid 4 He generated by ambient vibrations 
in the zero temperature limit (T < 100 mK) at various helium levels /iidc- The same data 
before (left) and after (right) image processing. Solid lines represent meniscus corrected 
eigenfrequencies in the annular volume according to Eq. CE3 with 7 = 375 pj/m 2 , L x = 
7r(25.0 mm + 27.6 mm)/2, L y = 1.3 mm, and ho = 1.75 mm. Periodicity of the mode marked 
with an ellipse is related to geometry 


the resonance width data shown in Fig. HOf ai can be somewhat extended at 
both ends when Af is too narrow or too wide to be determined directly. At¬ 
tenuation is discussed in details in another paper |Tj, where also the theoretical 
line shown in Fig. llOf al is explained. However, after revision described in the 
Ph.D. thesis by Matti S. Manninen, the theoretical estimate for the damping 
should have been smaller by factor of two, and the proper relation for the 
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/(Hz) 

Fig. 9 (Color online) Frequency spectra of the surface waves in superfluid 4 He excited 
by ambient vibrations at various temperatures. Helium level was h \dc =4.1 mm. 1 V 
corresponds to 0.17 mm. The color code represents the amplitude 



Fig. 10 (a) Resonance widths Af of the surface waves excited by ambient vibrations in 
superfluid 4 He at various temperatures at frequencies of / = 7.36 Hz (circles), 8.89 Hz 
(triangles) and 9.16 Hz (squares). Examples of the raw frequency spectra are shown in 
Fig. E] Solid line is the expected damping due to thermal quasiparticles. Helium level was 
/iidc =4.1 mm. (b) Amplitude dependency of Af with expected proportionality (gray 
lines). 1 V corresponds to 0.17 mm 
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resonance width is 

/ 15(pg + "/k 2 )h 3 u 4 ~ 15 pgh 3 u 4 

Q = Af = 8n 6 f(k B T) 4 ~ 8^f(k B T) 4 ' 

where u is speed of sound in 4 He. 


(17) 


6 Experiments on 3 He 

Experiments on superfluid 3 He are technically very challenging but bear in¬ 
terest due to the special properties of the superfluid states of 3 He. Our results 
here are much less systematic but interesting comparison with the results on 
4 He can still be made. 

It is not so easy to give proper readings for the temperature and the zero 
temperature limit is barely within reach at the zero pressure of our experi¬ 
ments on 3 He. We remind, that our cell was originally designed for studies of 
helium crystals, so that the installed quartz tuning fork thermometers were 
useless as they were located above the free surface of helium. The platinum 
NMR thermometer was only moderately thermally coupled to the experimen¬ 
tal volume and its readings saturated to a value around 0.4 mK. While the 
sample was cooled down by lowering the magnetic field B of the copper nu¬ 
clear cooling stage, the temperature was deduced from the demagnetization 
field assuming adiabaticity of the process, T ns oc B. 

The best independent reference point for such low temperatures is the 
superfluid transition temperature T c of 3 He, which could be observed by the 
interdigital capacitors themselves. When DC bias voltage was applied over the 
capacitor, the increase of the capacitance indicates that the higher the DC 
bias voltage the larger proportion of the surface of the capacitor was covered 
with helium both below and above T c . This effect was greatly enhanced in 
superfluid thus clearly showing the onset temperature. As a curiosity we point 
out, that even in an almost empty cell the IDC attracted substantial amount 
of helium remaining in the porous heat exchangers, when sufficiently high DC 
voltage was applied. 

The superfluid transition is demonstrated in Fig.QTla), where DC bias 
voltage was sinusoidally varied between —20 V and +20 V at 5 mHz. The 
measurement circuit with a 1 kHz capacitance bridge is shown in Fig. lllf b). 
The ferrite ring protected the capacitor against the noisy sensing channel of 
the bridge, which caused significant heating at the lowest temperatures with¬ 
out such filtering. Our value for the transition temperature T c = 1.07 mK is 
about 10% higher than often found from literature m- The transition tem¬ 
perature was deduced from the demagnetization field which was calibrated 
with the platinum NMR and other thermometers ( 3 He vapor pressure, cal¬ 
ibrated germanium resistor and noise thermometer) at higher temperature. 
Also our measured values for the phase transition temperatures T c , Tab and 
Tat at the melting pressure were about 10% higher than suggested by PLTS- 
2000 [12]. Such difference persisted despite of our most careful attempts to 
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perform proper calibrations for our thermometers using the accepted stan¬ 
dards, namely the vapor pressure of 3 He, as reference. 

Surface waves in 3 He were studied with three helium levels only since con¬ 
tinuous filling through the filling line at temperatures far below 1 rnK caused 
excessive heat load to superfluid 3 He and to the copper nuclear cooling stage. 
In order to fill the cell continuously, the compressible bellows volume should 
have been located below the experimental volume. 

Frequency spectra with the three 3 He levels are shown in Fig. fl2l together 
with calculated frequencies corresponding to the modes in the cuboid volume, 
Eq. (11611 . and the annular volume, Eq. m- The waves were excited either 
by ambient vibrational noise or by actively swinging the whole cryostat. As 
the resonances were measured only with three helium levels in rather complex 
geometry, it is difficult to reliably associate the measured and calculated fre¬ 
quencies with each other. The most distinct resonance at 16 Hz corresponds to 
m x = 3 in the cuboid volume though it was clearly seen only with the lowest 
helium level. 

As shown in Fig. [13] the surface waves in superfluid 3 He excited by ambi¬ 
ent vibrations of the environment vanished at temperatures T > T c /4 being 
observable again only above T ~ 50 mK. The temperature was deduced from 
the demagnetization field as the quartz tuning fork thermometers were located 
above the free surface. With full cell the lowest temperatures measured with 
the tuning forks were about 0.12T C at the melting pressure m ■ However, dur¬ 
ing the experiments on the free helium surface, heat leak to the experimental 
cell was reduced as the tube connecting the experimental cell with the bellows 
volume was empty. 




Fig. 11 (a) Temperature T ns deduced either from the platinum nuclear magnetic resonance 
thermometer (gray line) or from the magnetic field of the copper nuclear stage (black solid 
line). Superfluid transition temperature T c of 3 He serves as a fixed calibration point (dotted 
lines), which can be seen from the IDC data (oscillatory line) as helium is attracted towards 
high voltage capacitors more effectively below T c . (b) Circuit diagram to find T c 
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Fig. 12 Surface wave measurements in superfluid 3 He with three helium levels, hioc = 2.8, 
3.3 and 4.6 mm. FFT spectra generated by ambient vibrations (black lines). Fitted resonance 
widths and relative amplitudes of the actively excited waves (gray bars), see an example 
of the fits to raw data in Fig. Il5t a~). Thick gray lines represent eigenfrequencies in the 
cuboid volume without meniscus correction according to Eq. ©• Thick gray dashed lines 
are based on 2D FEM calculations according to Eq. (1161) with L x = 10.0 mm, 6 = 2.75°, 
and ho = 1.3 mm. Gray solid lines represent meniscus corrected eigenfrequencies in the 
annular volume according to Eq. tn with L x = 7r(25.0 mm + 27.6 mm)/2, L y = 1.3 mm, 
and ho = 1.75 mm. Here 7 = 155.7 pJ/m 2 
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Fig. 13 (Color online) Frequency spectra of the free surface waves in superfluid 3 He at 
various temperatures T ns of the nuclear stage. Here the waves were excited mechanically 
by noisy environment. Helium level was /iidc — 2.8 mm. 1 V corresponds to 0.2 mm. The 
amplitude is highlighted with the color code 


The resonance widths in superfluid 3 He at various temperatures are shown 
in Fig. ICTa'I together with the expected quasiparticle damping [1] 

n _ f _ (pg + 7k 2 )h 3 exp[E A /(k B T)\ pgh 3 exp[E A /(k B T)\ 
y Af 87 T*f P % ~ 87 T*fp% ’ 1 ' 

where E A and pf are the energy gap and Fermi momentum, respectively. The 
expected inverse proportionality between Af and amplitude is confirmed in 
Fig. Ildl bl . In contrast to 4 He, Af was wider than our frequency resolution 
even at the lowest temperatures. Thus there was no need to determine Af 
indirectly by integrating the FFT spectrum over the resonance. 

An alternative technique, phase sensitive active mechanical drive, can be 
used to study the temperature dependence of Af , see Fig. [15] The waves were 
excited with non-sinusoidal mechanical swinging and measured with higher 
harmonic frequencies. This method is very slow, since each individual fre¬ 
quency must be measured separately, see Fig. ESa). Temperature dependence 
of Af is weak [Fig. ITSl bl]. though the amplitude is temperature dependent 
even below 0.1T C [Fig. fTKI dl], This can be seen also in Fig. [T5T cl where Af 
saturates from the expected inverse proportionality to the amplitude at the 
lowest temperatures. In addition, the frequency widths of the modes with dif¬ 
ferent frequencies saturate at different temperatures to value about 0.1 Hz 
for most of the modes. In principle, the heat load caused by the mechanical 
excitation could depend on frequency, though with full cell the superfluid 3 He 
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Fig. 14 (a) Resonance widths Af of the surface waves excited by ambient vibrations in 
superfluid 3 He at various temperatures T ns at frequencies of / = 5.04 Hz (circles), 5.41 Hz 
(triangles), 6.50 Hz (squares) and 8.57 Hz (diamonds). The raw spectra are shown in Fig. |l3l 
Helium level was /ijdc = 2.8 mm. (b) Amplitude dependency of Af with expected propor¬ 
tionality (gray lines). 1 V corresponds to 0.2 mm 


temperature measured directly with quartz tuning fork thermometers was not 
excitation frequency dependent. In Fig. [T5T df the wave amplitude is scaled by 
frequency, which corresponds to the vertical velocity of the wave. Even at the 
lowest temperatures the wave velocity is very low, vpf -C ksT, which was 
assumed in Ref [1) when calculating the damping due to quasiparticles shown 
in Fig. fl5lbl. 

For comparison with the data measured in superfluid 3 He, Fig. [T6] illus¬ 
trates the measurements in 3 He above the superfluid transition temperature. 
Here again the waves were excited by swinging the whole cryostat pneumat¬ 
ically. The temperature dependence is seen merely in amplitude [Fig. flTTlf cl] 
while Af is either temperature independent or the dependence is weak [Fig. 
Mb)}- This is not what was expected. As for the waves excited by ambient 
vibrational noise, the prediction was, that Af and amplitude would have had 
inverse correspondence. 


7 Conclusions 

Surface waves were studied both in superfluid 4 He and 3 He. Superfluid 4 He 
level was raised continuously below 100 mK and 12 frequency modes up to 
60 Hz were observed corresponding to the geometry of the main cuboid ex¬ 
perimental volume. In addition, there were several modes originating from the 
surrounding annular volume coupled to the cuboid volume via two narrow 
channels. In superfluid 3 He the highest resonance frequency was 16 Hz cor¬ 
responding to the third mode in the cuboid volume. The resonances in the 
annular volume could not be identified as the waves were studied only with 
three distinct 3 He levels. Moreover, the characteristics of the frequency spec- 
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Fig. 15 (a) Phase-sensitive detection of surface waves in superfluid 3 He actively generated 
by swinging the cryostat. Detection frequency / was five times higher than the excitation 
frequency. Black solid lines are fits to the measured data (black dots). Resonance widths 
and amplitudes of the fitted curves are visualized as gray bars. Here Hidq =4.6 mm and 
the temperature of the nuclear stage was in the zero temperature limit for the superfluid, 
Tns = 0.03T C . (b) Temperature dependency of resonance width Af. Solid line is the expected 
damping due to thermal quasiparticles, (c) Amplitude dependency of Af. (d) Temperature 
dependency of the wave amplitude a at a resonance. The amplitude a is multiplied by 
resonance frequency / which corresponds to vertical surface velocity on the IDC 
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Fig. 16 (a) Phase-sensitive detection of surface waves in normal 3 He actively generated 
by swinging the cryostat. Detection frequency / was three times higher than the excitation 
frequency. Black solid lines are fits to the measured data (black dots). Resonance widths 
and amplitudes of the fitted curves are visualized as gray bars. Here T = 75 mK and /iidc = 
3.1 mm. (b) Temperature dependency of resonance width Af. (c) Vertical RMS speed of 
the resonance on the IDC deduced from the RMS amplitude a multiplied by resonance 
frequency / 


tra in 4 He and in 3 He differed significantly though they were supposed to be 
similar oscillations in the same geometry. 

The surface wave resonance frequencies were measured with great preci¬ 
sion especially in superfluid 4 He. In principle surface tension could have been 
deduced from those frequencies. However, this is not straightforward since, 
besides our nontrivial cell geometry, the effect of surface tension is twofold. 
In addition to the direct surface potential energy it has an additional effect 
to the surface shape by forming menisci near the walls. The measurements 
were in reasonable agreement with a 2D numerical model where both effects 
were taken into account. Also the basic analytical model fits to the data when 
kinetic energy is scaled by factor of (1 + 0.014m x ) to effectively account for the 
meniscus correction. Better agreement is expected with a full 3D numerical 
model where surface tension and porosity of the sintered layers are properly 
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taken into account. This may also result a more precise value for the surface 
tension as the discrepancy among previous measurements is about 6% mm 

The temperature dependence of the resonance width was studied with two 
different methods. When the waves were generated by ambient vibrational 
noise, the resonance width was inversely proportional to amplitude as ex¬ 
pected. However, for the waves in 3 He actively generated by swinging the whole 
cryostat the temperature dependency of the width was weak even though the 
amplitude was temperature dependent. Above T c the width was essentially 
temperature independent. 

Admittedly, our experimental arrangement was not optimal for producing 
easily interpretable data on free surface resonances in superfluids, as it was 
initially designed for studies of melting-freezing waves, or crystallization waves, 
on the superfluid-solid interface. Since the experiments performed at very low 
temperatures are usually not easily reconfigured, any modifications to the 
present system were not undertaken. Further experiments using a more suited 
setup designed particularly for this purpose would probably clarify some of 
our unfit results. The presented study was performed, nevertheless, since it 
was a straightforward extension to the solid phase experiment and since very 
few or no similar studies exist. Moreover, the basic features are similar for the 
waves both on the free surface and on the mobile superfluid-solid interface. 
Further discussion on both phenomena can be found from the Ph.D. thesis by 
Matti S. Manninen. 
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